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1 Introduction 

^ Let K C M"^, (i > 2, be a compact convex set with interior points and let be a Poisson 

I point process in M with intensity measure /i = Xid with A > and the (i-dimensional 

j> Lebesgue measure £d- For every point x E r] we define the Voronoi cell of x by 



X 



= [z e R'^ : \\x - z\\ < \\y - z\\ for all y e r]] 

and call x the nucleus of V^. We have int(\4) fl int(Vy) = for x 7^ y G r] and 
[Jx&ri = R'^ such that the collection (V^)^g^ of random polytopes constitutes a ran- 
dom tessellation of M*^, the so called Poisson-Voronoi tessellation, which is one of the 
standard models in stochastic geometry, and we refer to [12] and the references therein 
for further details. 

For our set K we define the Poisson-Voronoi approximation A{K) as 

A{K)= U Vx, 

xGrinK 
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Figure 1: Poisson-Voronoi approximation of an ellipse 



which is a random approximation of K. It is possible to interpret the Poisson-Voronoi 
approximation in the following way. One wants to reconstruct an unknown compact 
convex set K, but the only information available is a kind of oracle which tells for 
every point of a realization of the Poisson point process if it belongs to K . Now one 
approximates the unknown set K by taking the union of the Voronoi cells with nuclei 
in K. 

In this paper, we are interested in the volume of the Poisson-Voronoi approximation 

PV(K) = Vol(A(ir)). 

A short computation yields EPV(-ft') = Vo\{K), which means that PV(-ft') is an un- 
biased estimator for the volume of K. Under weaker assumptions than convexity 
on the approximated set K, it was shown in [U El |H] that PV(fr) — Yo\{K) and 
Yo\{A{K)^K) — > as A — 7- oo, where A{K)AK stands for the symmetric difference 
of A{K) and K. In [2], upper bounds for the asymptotic behavior of VarPV(ii') and 
Var Vol(A(ii')Aii') and large deviation inequalities for PV(iir) and Vol(A(i^)Ai^) are 
derived for the same setting as in this paper. In [10], these results are extended to all 
non-centered moments and to a more general class of approximated sets, namely sets 
of finite perimeter. 

The main result of this paper is that 'PY{K) behaves asymptotically like a Gaussian 
random variable if the intensity of the Poisson point process goes to infinity. 

Theorem 1.1 Let N he a standard Gaussian random variable. Then 
PV{K) - Yo\{K) 



v/Var PV(A"; 



N in distribution as A — > oo. 



As it is pointed out in [2j, the Poisson-Voronoi approximation has applications in 
nonparametric statistics, image analysis and quantization problems. In this context. 
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Theorem can be helpful, since it allows to treat PV(i^) as a Gaussian random 
variable if the intensity of the Poisson point process is sufficiently high. 



For the proof of Theorem 1.1 we use so called Wiener-Ito chaos expansions, which 
give us a representation 

oo 

py{K) = yo\{K) + Y,Wn), 

n=l 

where the functions /„ G -/^^((M*^)"), n G N, are known and /«(■) denotes the n-th 



multiple Wiener-Ito integral. The key argument of our proof of Theorem 1.1 



IS an 



abstract central limit theorem, which is derived in Section M and could be helpful for 



other problems as well. This Theorem 3.1 is based on a general result for the normal 



approximation of Poisson functionals (see Theorem 3.1 in [B]), which is used in a similar 
way as in p]. 

In order to check the assumptions of this abstract central limit theorem, we have 
to prove some kind of uniform convergence for n!||/„||^. Combining this property with 
the identity 

oo 

VarPV(i^) = 5^n!||/„||^, 

n=l 

we obtain clS db byproduct: 

Theorem 1.2 Let tk he the inradius of K. Then, there are explicit constants C_,C>0 
depending only on the dimension d such that 

d-l 

Ck,V,.^{K)X-'--^ < VarPV(ir) < CJ2^'d-^V^{K)X-^+'^ (1) 

for A > {2/rxY, where Vi{K)^i = 0,...,n, are the intrinsic volumes of K and kj 
stands for the volume of the unit ball in W . 

It is well known that VaiK) = Vo\{K), Va-iiK) = lS{K), where S{K) is the surface 
area of K, and Vo{K) = 1. Both bounds in ([l| are of order A~^~3 such that Var PV(i^') 
also has order A^^^^. 

The upper bound in ([T]) is also contained in [2], where it is proven by a combination 
of the theory of valuations and the Poincare inequality. The Poincare inequality is also 
connected to Wiener-Ito chaos expansions (for more details we refer to [5j). The lower 
bound is new as far as we know. 

Although the construction of the Poisson- Voronoi approximation does not depend 
on the convexity of K and can also be done for more general classes of sets, we formulate 
our main results only for convex sets, in order to simplify the proofs. In the final 
Remark |4| we give two alternative assumptions for the approximated set that allow us 
to replace the convexity assumption. 

This paper is organized in the following way. In Section [2| we introduce Wiener- 
Ito chaos expansions and recall the coarea formula we use several times in our proofs. 
An abstract central limit theorem for Poisson functionals and a helpful proposition to 
check one of the conditions are derived in Section [3j In Section |4], the Wiener-Ito chaos 
expansion for the volume of the Poisson- Voronoi approximation is computed and used 
to establish Theorem L2 before the proof of Theorem 1A_ is concluded in Section [5] 
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2 Preliminaries 



Wiener-Ito chaos expansions. Our main tool in this paper are so caUed Wiener- Ito 
chaos expansions, which are briefly introduced in the following, and we refer to [HI EJ [7] 
for more details. For a Poisson functional F = F{ri) depending on a Poisson point 
process rj over a Borel space (X, A",//) with a cr-flnite non-atomic intensity measure /z 
(for PV(i^): X = M'^, X is the standard Borel cr-fleld in M*^ and (x = \id) one deflnes 
the difference operator as 

D^F = Fir] + 6,) - F{ri) 

for X G X, where 6x is the Dirac measure concentrated at the point x. The difference 
operator has the geometric interpretation that it measures the effect of adding the 
point X to the Poisson point process rj. Therefore, it is sometimes called add-one-cost 
operator. The iterated difference operator is given recursively by 

■D Xl,...,XnP F)xiF)x2,...,Xn-F (^) 

and is symmetric under all permutations of Xi, . . . , Xn- The deflnition in ^ is equiva- 
lent to 

Dx„...,x„F= Yl (-irlV(r^ + 5^5.J. (3) 

IC{l,...,n} i&I 

Denoting the n-th multiple Wiener-Ito integral with respect to the compensated Pois- 
son point process 77 — /i by /„(■) and deflning functions : X" — )■ M := M U {±C)o} 
as ^ 

fn{Xi, ...,Xn) = —EDx„...,x„F, (4) 

we have the following representation for square integrable Poisson functionals F (see 
Theorem 1.1 and Theorem 1.3 in ^): 

Theorem 2.1 Let F G ^^(P). Then /„ G ^^(X") forneN and 

00 

F = EF + 5^ /„(/„). (5) 

n=l 

Moreover, 

00 

VarF = ^n!||A||2, (6) 

n=l 

where \ \ ■ ||„ stands for the usual norm in L^(X"). 

We call the identity ^ the Wiener-Ito chaos expansion of F and the functions fn the 
kernels of the Wiener-Ito chaos expansion of F. Note that ([s]), which is sometimes 
denoted as Fock space representation, is an orthogonal series since 



EI^IM 



n\ fg (i/i", n = m 
0, n ^ m 
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Analogously to Q the covariance of two Poisson functionals F,G G L^(P) with F = 
+ En=i Wn) and G = EG + ^^=1 UQu) is given by 

oo 

Cov(F,G) = 5^n!(/„,^„)„, (7) 

n=l 

where (■, ■)n is the usual inner product in L^(X"'). 

In the next section, we need the inverse Ornstein-Uhlenbeck generator L^^ which 
is for centred random variables with a chaos expansion (|5| given by 

oo ^ 

L-'F = -Y,-Wn). (8) 

n=l 

In this context, it is also possible to define the difference operator as 

oo 

D,F = ^n/„_i(/„(x,-)) 

if Er=i^^!|l/n||^<oo. 



n=l 



Coarea formula. Our main tool for the computation of integrals where the kernels 
of the chaos expansion of F arise, e.g. ||/n||^, is the so called coarea formula. By "H*", 
we denote the m-dimensional Hausdorff measure. If / : — )> M" is differentiable in 
X G M™, we define the Jacobian J f\x) by 

J/(x) = v/det(/'(x) r{xY), 

where /' stands for the Jacobi matrix of /. Note that a Lipschitz function is almost 
everywhere differentiable such that its Jacobian is almost everywhere defined. Using 
this notation, we have (see Corollary 5.2.6 in P], for example): 

Theorem 2.2 // / : — )■ M" a Lipschitz function and m > n, then 

[ g{x)Jf{x)Udx) = [ [ g{z)n"'-^{dz)iM) 

holds for each Lebesgue measurable B C M™" and each nonnegative irn-n^'^O'Surable func- 
tion g : B ^R. 

For n = 1, we have Jf{x) = ||V/(a:)||. Note that || ■ || stands for the usual Euclidean 
norm, whereas || ■ is the norm in L^(X") or L^((]R°')"'). 



3 An abstract central limit theorem 

In this section, we prove an abstract central limit theorem for a more general setting. 



which is used to show Theorem 1.1 As in the previous section, we assume that r/ is a 
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Poisson point process over a Borel space {X, X, fi) with a a-finite non-atomic intensity 
measure and F e L^(P) is a Poisson functional with a Wiener-Ito chaos expansion 

oo 

F = EF + J2Wn). 

n=l 

We are interested in the Wasserstein distance between our Poisson functional and a 
Gaussian random variable, which is very helpful for establishing central limit theorems 
since convergence in Wasserstein distance implies convergence in distribution. The 
Wasserstein distance dw of two random variables Y and Z is given by 

dw{Y,Z)= sup \Eh{Y) -Eh{Z)\, (9) 
/leLip(i) 

where Lip(l) is the set of all functions /i : M — > M with a Lipschitz constant less or 
equal than one. In this setting, we can state the following: 

Theorem 3.1 Let F G L^(P) and let N be a standard Gaussian random variable, 
a) For every A; G N one has 



, F-EF \ vE„=fc+i«!||/n|L , ^ JE~j ,J^^. 

dw , =,iV < 2-^ ; + fc > ^^^ + fc2> 10 

VVVai^ v/Vai^ VarF ^VarF ^ ^ 

i?., = E(/,_i(/,(z, ■)), h^x{JA^. ■))? - (E(/.-i(/.(z, ■)), ■))))' 
= E(/,_l(/i(^,■))^/.-l(/.(^,•))') 

jor i,j = l,...,k. Here (-, ■) stands for the usual inner product in L^(X). 
b) Let F depend on a parameter A > 0. // there are constants c„ G M such that 

n\\\f IP 

^^^p < On for all A > Ao G M and ^ c„ < oo (11) 

and 

^fR~, as A^oo (12) 



VarF VarF 
for all i,jEN, then 



,'F-EF 
dw , — )■ as A — 7- cxD. 

VVa^F 



Since F is standardized in Theorem |3.1[ we can assume without loss of generality that 
EF = 0. The idea of the proof of ( |10[ ) is to define truncated Poisson functionals F^, 

G N, by 



n=l 
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and to use the triangle inequality 
F 



d 



w 



VVaFF' 



N \ <d 



w 



WaFF' WaFF 



+ d 



w 



N . 



(13) 



Now we compute upper bounds for both expressions on the right hand side in (13). 
Lemma 3.2 It holds 

F F, 



d 



w 



WaFF' WaFF 



< 



X]n=A:+l ^-1 1 /" I Ira 

VarF 



(14) 



for every A; G N. 



Proof. By the definition of the Wasserstein distance in ^ and the Cauchy-Schwarz 
inequality, we obtain 



d 



w 



WaFF' VVaFF 



sup 

/leLip(i) 



VVaFF 



VVarF VVarF 



v^VaFF 
E|F - Ffel 
/VaFF 



< E 



^ v/E(F-F,)^ ^ VE(E:r=.+iWn))' 



WaFF 



VVaFF 



and ( 14 ) is a direct consequence of (|6| . 



□ 



For the second expression in ( 13 ) we need the following inequality: 



Lemma 3.3 We have 



d 



w 



WarF 
for every A; G N. 



X]n=fc+1 1 1/™ I In 



VarF 



l<i,j<fc 



VarF 



i=l 



VarF 



(15) 



Proof. Theorem 3.1 in [6] tells us that 
Fk 



d 



w 



v^VaFF' 



< E 



+ 



1 - ^ / D^Fk i-D^L-'Fk) ^{dz) 
E[{D,Fkf\D,L-'Fk\]fi{d. 



(VarF) 2 Jx 



where L ^ is the inverse of the Ornstein-Uhlenbeck generator as given in dsj). By 



the fact that 1 = Z]n=i "-'ll/nlln/ VarF + ^^^^^ n!||/„||„/ VarF and the triangle 
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inequality, we obtain 



d 



w 



V VVaiF 



N < 



Sn=fc+1 I In 



VarF 



+ 



1 



VarF 
1 



y;n!||/„||^^- / D,Fk{-D,L-'Fk)fi{d2 
E[{D,Fkf\D,L-^Fk\\^l{dz). 



(VarF)2 Jx 

In the proof of Theorem 4.1 in |9j, it is shown that the last two expressions are bounded 
by 



l<i,j<k 



i=l 



VarF' 



which leads to (15). 



□ 



Proof of Theorem 3.1, The inequality (10) in a) is a direct consequence of (13), 
Lemma 3.2 and Lemma 3.3 If the conditions of b) are satisfied, for every e > we 
can find ko such that 



Sn=fco+l '^'ll/nlln 



VVarF 



< 



for all A > Aq. Because of (12) it exists a constant A > such that 

Ri 



VarF 



i=l 



VarF 



for all A > A. Combining these inequalities with (10), we obtain 

'F-EF 



d 



w 



VVar F 



for all A > max {Aq, A}. 



□ 



Remark 1 Our abstract central limit Theorem 13.11 and Theorem 11.11 have the draw- 
back that they do not give a rate of convergence. This problem is caused by the 
truncation step. The second expression in (13) vanishes as A — oo. But the first 
summand does not necessarily converge for a fixed as A — )■ oo. By taking — )■ oo in 
the previous proof, we obtain convergence as A — )■ oo, but cannot give a rate of conver- 
gence. An alternative approach would be to apply the underlying general central limit 
theorem (see Theorem 3.1 in ^) directly to the infinite Wiener-Ito chaos expansion, 
which gives a sum of an infinite number of expected values of products of multiple 
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Wiener-Ito integrals as an upper bound. These summands are also the Rij and Ri in 



our Theorem 3A and it is possible to show an upper bound and a rate of convergence 
for each of them as one will see in Section |5} But it seems very hard to prove that the 
series over these bounds converges. 



In order to neatly formulate a very helpful criterion for the condition (12) in The- 
orem 3.1, we need the following notation: Let Ilj , be the set of all partitions tt of the 



variables 

^(1) ^(1) ^(2) J2) (3) (3) (4) (4) 

•^1 J ■ ■ ■ J -^i ) -^1 ) • ■ • ) -^i ? -^1 ) • • • ) j ) -^1 T ■ ■ ■ 1 -^j 

such that 

• all variables with the same upper index are in different elements of vr 

• every element of vr has at least two variables as elements. 

By Hi J, we denote the set of all partitions vr in Uij such that for any decomposition of 
{1,2,3,4} in two disjoint nonempty sets Mi and M2 there exist h G Mi and I2 G M2 
such that two variables x^^^^ and x^^f are in the same element of vr. 

We say that a partition vr connects two variables if they are in the same element of 
71, and |7r| stands for the number of elements of a partition vr G Ilij. By {fi* fi* fj* fj)-^-, 
we denote the function from X''^' to M we obtain if we replace all variables that are in 
the same element of tt by a new variable. Using this notation, we can give the following 
upper bounds for R^j and Ri. 

Proposition 3.4 // \ {fi * f^* f^* fi)^\d^\'^\ < 00 for all tt G 11^,^ and 
IxM K/i * fj * fj * /jO'tM/^'''' < 00 for all n G Ujj, it holds 

R^J< T I \Ui*h*f3*f3)M-^- (16) 

Moreover, one has 

R^<y2^ \{f^*f^*f^*fi)W''^. (17) 
TreU-i.i 



Proof The fact that the integrals over (/j * fi * fi * fi)n and {fj * fj * fj * fj).,, are 
finite ensures that we can apply the product formula for multiple Wiener-Ito integrals 
(see Theorem 3.1 in fH] or Proposition 6.5.1 in [7]) to Jj_i(/i(s, •))/j_i(/i(t, ■)) and 
Ij-i{fj{s, ■))Ij-i{fj{t, ■)) for /i-almost all (s,t) G X^. This formula gives us the ker- 
nels of the chaos expansions of Ii_i{fi{s, ■))/i_i(/j(t, ■)) and Ij-i{fj{s, ■))/j_i(/,(t, ■)). 
Combining this with the covariance formula ([T]), we know that 

E / /" /,_i(/,(s,-))/.-i(/.(t,-))Vi(/.(^'-))Vi(/.(^'-))M^%(^^) 
Jx J X 

= nil! (h{s,-)* f,{t,-)* f,{sr)* f,{t,■))^d^^\-\^^{ds)^^{dt). 
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For i ^ j there is no partition connecting either variables of the first and the third 
function or variables of the second and the fourth function. Hence, we have only 
partitions from 11, j if we add s and t to the partitions. For i = j we have — 
partitions where the variables of the first and third and of the second and fourth 
function are pairwise connected. The sum over this partitions is {i — — l)!||/j||^. 
But exactly this term is subtracted for i = j. If we add s and t to the remaining 
partitions, we also obtain partitions from Ilj^j. For Ri we can also apply the product 
formula and it holds 



X 



{fi{z, ■) * fi{z, ■) * fi{z, ■) * fi{z, ■))^dii\'^\ii{dz) 




X J x\^\ 



^ JxJxM 



□ 



4 Proof of Theorem 1.2 



Because of Theorem 1 in [2], we know that PV{K) G L^(P) and Theorem 2.1 implies 
the existence of a Wiener-Ito chaos expansion. In the following, we compute the kernels 
of this decomposition and use ^ to prove our bounds for the variance of PV{K) in 
Theorem 11.21 



From now on, we denote by p(-, ■) the usual Euclidean distance, which is given by 
p(x, y) = \ \x — y\ \ for two points x, ?/ G M*^ and p(x, A) = inf p(x, y) = inf I |x — I for 

yeA y&A 

X and A C W^. Moreover, B'^{z, 6) stands for a ball with center z and radius 6 in 
M''. Using this notation, we have the following formula for the kernels of the Wiener-Ito 
chaos expansion of PV{K): 

Lemma 4.1 Letxi, . . . ,a;„ G M*^. Fory G we define x{y) := aig'ma.Xx=xi,...,xn piv, x) 
and z{y, rj) := arg min^g^ p{y, z) . Then 

= U{J^^y)(^K)nz{y.v)iK''^B\y^\y~x{y)\\))dy {l^^ 



-1 



l{x{y) G K)nz{y. v)iKU B\y, \ \y - x{y) \\))dy. 



Proof. Since rf) is the nucleus of the Voronoi cell y belongs to, it is easily seen that 
PV(K) = Yo\{{y G : z{y, r]) e K}) = [ r/) G K)dy. 
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Combining this with ([3]), we obtain 

lG{l,...,n} iel 



[ V i-ir+\'\liziy,vl^{x.:teI})eK)dy. 



IC{l,...,n} 

Now we consider the sum of the indicator functions on the right hand side for a fixed 
y E K. Let imax be the index of the Xi that maximizes x{y). For J C {1, . . . , n} \ {imax} 
with / 7^ 0, it holds z{y,ri U {xi : i G /}) = z{y,ri U {xj : i G / U {imax}}) and the 
summands for / and I {imax} on the right hand side cancel out because of the different 
signs. Hence, we obtain 

D,,_,^FViK)= [ i-iniiziy,r])eK)-l{z{y,r]U{xiy)})eK))dy. 

Now it is easy to see that 

liziy,r]) eK)~ l{z{y,r]U{xiy)}) G K) 

1, p{y, x{y)) < p{y, z{y, r])), z{y, r]) G K, x{y) ^ K 
-1, Piy,x{y)) < p{y,z{y,r])), z{y,ri) ^ K, x{y) e K 

Combining this with the definition of the kernels in (|4]), we obtain 
Ux,,...,Xr,) = / l{x{y)^K)F{z{y,r^)^K^UB''{y,\\y-x{y)\\))dy 



-r 



nl 



l{x{y) e K)F{z{y,r]) ^ K U B\y,\\y ~x{y)\\))dy. 

□ 



Remark 2 For fi we have the representation 

'eVo1({z G : p{z,x) < p{z,r]nK^) < p{z,r]nK)}), xeK 
-EVol{{z eR"^ : p{z,x) <p{z,rinK) < p{z,r] n K^)}), x e 



/i(^) 



which means that |/i(x)| is exactly the volume of the points that change between A{K) 
and A{K)'-' if the point x is added to the Poisson point process. 

Our next goal is to compute upper bounds for ||/„||^ such that we obtain by ([6| an 
upper bound for the variance of PV(i^) and can check condition (11) in Theorem 3.1 In 
formula (18), the distance between a point y G M*^ and x{y) plays an important role. In 
order to handle this quantity in the following, we define functions hn : (M'^)" — t- M x M*^ 
by 

hn{xi, . . . ,Xn) = {min max p(?/, Xj), argmin max p{y,Xi)). 

yeR i=l,...,n ydiR i=l,...,n 

From a geometrical point of view, hn gives the radius and the center of the smallest 
ball that contains all points 

The function allows us to give the following upper bound for /„: 
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Lemma 4.2 Letxi,...,Xn G and let r = hn\xi, . . . ,Xn) = niin max p{y,Xi). 

j/SM i=l,...,n 

Then ^ 

\fn{xi,...,Xn)\ < Yyy^exp(-A/td/). (19) 



Proof. As a consequence of Lemma 4.1, one has 



|/n(xi,...,a;„)| < ^ / nziy,r])^B''{y,\\y-xiy)\\))dy. 

By the definition of r, we know that the sets M.'^ \ B'^{xi, r),i = 1, . . . ,n, cover 
Combining this with the previous inequahty and using polar coordinates, we have 

\fn{xi,...,Xn)\ < nz{y,ri)^B''{y,\\x,-y\\))dy 

^- JR<*\B<*(xi,r) 
I ^ r 

= —,y2 exp{-\Kd\\xi- y\\'^)dy 



n 

1 



(n- 1)!A 



exp(-A/trfr'^)r'^ ^dr 
exp(-Aftrfr'^). 



□ 

By definition, fn{xi, ■ ■ ■ , Xn) measures the effect on PV(-ft') of inserting points. Lemma 



4. 2| reflects the fact that this effect is small if the distances between the points are large. 
Similar one expects that /„(a;i, . . . , x„) is small if all points are close together but are 
far away from the boundary of K. This effect is described by the following Lemma: 

Lemma 4.3 Let Xi, . . . , x.„ G and {r,y) = hn{xi, . . . , x„). If 5 = p{y, dK) > 8r, 
then 

|/.(xi,...,x„)| < ^exp{-XKj''/8''). (20) 

Proof. Since p{y, dK) > 8r, all either in K or K^. Let x = ^{y + 

projgj^(y)), where projg^(^) stands for the metric projection of y on the boundary of 
K. If ^ G K, it can happen that the metric projection on dK is not unique. In this 
case, it does not matter which of the points is taken. Then, we have 

I < p{x, y) < ^5 < p{y, dK) for all y G B''{x,, 6/8) 
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and, by (18), it follows 



\fn{Xi, 



^- J R<*\Brf (11,5/8) 



[z{y,r^)iB\y,\\y-x,\\))dy 



+ - / ¥{z{y,r^)iB\y,p{y,dK)))dy 

< \ I n<y,v)iB\y,\\y-xMdy 
+ -, I ¥{z{y,ri)iB\y,\\i-y\\))dy. 

^- jR\B'i{x,6/8) 



A straightforward computation as in Lemma 4.2 yields (20) 



□ 



Combining Lemma |4.2| and Lemma |4.3| leads to the bound 

...,Xn)\< J^{hn{Xi, . . . , X„)), 

where : M x M'^ — M is given by 

7 (r = I (^^^ exp(-A/trfr'^), p{y, dK) < 8r 
X^^^vi-^^^ar'^), piy,dK)>8r 

By the coarea formula Theorem 2.2[ we obtain for n > 2 

ll/nlln < A" / J^{hn{Xi, . . . ,Xn)ydXi . . .dXn 



(21) 



= AM / / f,Xr,yfJhn{x,,...,Xn)-'n^'-''-\d{xi,...,Xn))dydr 

Jo JR^ Jh~^{r,y) 

POO P P 

= A" / / 7„(r,y)M J/i„(xi,...,x„)-iH"''-'^-^(d(xi,...,x„))rfyf/r. 

Jo Jr'' Jh~^(r,y) 

It is easy to see that hn{axi + v, . . . , axn + v) = a/i„(xi, . . . , x„) + (0, v) for all a > 
and V e M*^ and a short computation shows h'^{axi + t>, . . . , aXn + v) = h'^{xi, . . . , Xn), 
which implies 

Jhn{aXi + V, . . . , aXn + V) = Jhn{Xi, . . . , Xn) 

for all a > and f G M'' and 

[ JhniXi, Xnr'W"'^''-\d{Xi, Xn)) 

JK{x,, x„)-^H"'^-"-i(d(xi, . . . , Xn)). 



r 



y) 

nd—d—l 



Hence, (21) simplifies to 

1 1 /n I In ^ C'nA" 




f^{r,y)it,y)--'-'/''dydr 



(22) 
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for n > 2 with constants 



/ln'(0,l) 



Lemma 4.4 The constants Cn, n > 2, are finite and there is a constant q > only 
depending on the dimension d such that 



n 



a<c,( . )(„-l) 



for n> d + 1. 



Proof. A straightforward computation yields 

Cn — ^ Jhn{Xi, . . . , Xn) "H {d{xi, . . . , Xn)) 

Jh-\o,i) 

= (n - / / / Jhn{xi,...,Xn)-'n'"'~''-\d{xi,...,Xn))dydr 

Jo JB^iO,!) Jh-\y,r) 

= {n - l)dK~''^^^'^ l{hnixi,...,Xn) e[0,l]x B'^{0,l))dxi...dxn<oo. 

J{R'i)" 

For almost all {xi, . . . , Xn) € (M'^)" at most d + 1 points are on the boundary of the 
minimal ball that contains all points, and we assume that these are Xi, . . . , x^+i- Since 
the center of the minimal ball is in B'^{0, 1) and the radius is in [0, 1], these point must 
be in B'^{0, 2). The remaining points are in a ball with radius 1 around a center given 
by the first d + 1 points. These considerations lead to the bound 



n 



for n > + 1 . 



□ 



Our main tool for the computation of the right hand side in (22) is the following 
inequality: 

Lemma 4.5 There are constants Ci^^, C2^d > only depending on the dimension d such 
that 



d-l 



fn{r,yfdy < exp{-2XKy)J2^^d-^Vi{K) 
for all r > and n >2. 



1=0 



CiAiKgr'^y d C2,d{l^dr'^) 
((n-l)!)2A2 (n!)2A3 



Proof Let {dK)s = {y eR'^ : p{dK,y) < s}. Together with ^ and we obtain 
7nir,y?dy < [ \—:^expi-2XK,r'')dy (23) 



J(dK)sr {{n-iy.yx^ 

+ [ —^expi-2Xf,,p{y,dKy/8'')dy. 
14 



The volume of {dK)sr H K'-' is given by the Steiner formula (see Theorem 2.2.4 in [TT] ) 
and Vol{{dK)sr n K) < Vo\{{dK)8r n K^) such that 

J{dKw [{n-iy.yy 

^ d-l 



< 



i=0 
d-l 



exp(-2A/€,r'^) 5^ /«,_,y,(ir)(/s:,r'^)i-i 



((n-l)!)2A2 



To the second expression in (23) we apply the coarea formula with the Lipschitz- 
function m : M'^ — )■ M, x t-)- p{x, dK). 

^ exp{-2XKdp{y,dK)^/8'')dy 




^ exp(-2A/«d578'^) \\Vu{y) \\-^n'^'\dy)d6. 



8r 



(n!)2A2 



It is easy to see that |V^(2^^)M(a;)| = ||i;(x)||, where Vv(x)u{x) is the directional deriva- 
tive in direction v{x) = x — proj^'(a;). Hence, we have ||f(a;)|| = |V„(2;)M(a;)| < 
||V'u(x)|| ||f(a;)|| and ||Vu(x)|| > 1. By the Steiner formula (see Theorem 2.2.4 in 
ITT]'), we know that 



n''-\{z G R'' : piz, dK) = 6}) < 2H''-\{z E K^' : piz, dK) = 5}) 

d-l 

= Y,(d-i)t,,.,V,{K)5''--^-' 

i=0 

and altogether we obtain 

/ exp(-2A/trfp(y, 9A078")f/2/ 

/CO A 
— — exp(-2A«:,(578'^) V(rf - i)Kd-^V,iK)5''-'-'d5 
.r (n!)2A2 ^ 

^ n\^E(^-^)'^<^-^*W(80~W exp(-2A«:,578'^)5'^-M5 

od d.-^ 

= 2n\^ exp(-2A«,r'^) V(d - z)/€,_,F,(ir)(8r)^* 



< 



C24 



i=0 
d-l 



(n!)2A3 



exp{-2XKy)J2^^d-^V.{K){K,r'')-^^. 



1=0 



□ 

By hi{xi) = (0,xi), Lemma 4.3 and the coarea formula with the same function u as 
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in the previous proof, it follows that 

ll/ill? < A / ^exp(-2A«:,p(i/,9ir)V8'^)cil/ 

exp(-2A/€,r78^) \\Vu{y) \\-^U''-\dy)dr 



(24) 



< 



4 

A Jo J«-i(r) 
4 



4 ;"00 

- / exp(-2A«dr78^) - ^)/€,_,V,(ir)/-i-^rfr. 



Combining (22) and Lemma 4.5, we have 

d-l 



ll/nll^ < Cnl exp(-2AKrf/) VK,-,V-,(ir) ^)'" ' ( A/s:,r")"- ^ 



«QQ d 1 

+a / exp(-2A/€rfr'^) V K,^.y^{K) 

■^0 i=o 



C2,dA 



-2+5 



-(A/trfr 



dNn— 1- 



<* dr. 



for n > 2. Comparing this with (24), we see that the same bound also holds for n = 1 
if the constant Ci is chosen appropriately. Now substitution and the definition of the 
Gamma function lead to 



d-l 



A-2+^ 1 



i=0 



2"- d (2AKrf)^ Jo 



exp(-?/)?/" ddy 



d-l 



A-2+^ 1 



i=0 



d {2XKd)d Jo 



exp(-y)y" ^ ddy 



d-l 



i=0 



Ci^Tin - ^) ^ C24T(n - 1 - 



((n- 1)!)22" 



+ 



(n!)22" 



A-2+3 



with constants Ci^^, C24 > and it is easy to see that 



\ ( ^ . 1 A "'""^ 

nil' < (^5Mn + 52,d ^^^_^J 5^/€rf_,F,(ir)A 



(25) 



By Lemma |4.4[ we know that C„ is bounded by a polynomial of order d + 2 in n and 
it follows directly that the series 



n=l ^ 



' n{n — 1) 



converges, which proves the upper bound in Theorem 1.2 and that the condition (11) 
in Theorem 3A is satisfied for the Poisson-Voronoi approximation. 

In order to conclude the proof of Theorem L2, it remains to construct a lower 
bound. Because of ^ it is sufficient to give a lower bound for ||/i||i. 
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Lemma 4.6 There is a constant C_ only depending on the dimension d such that 

\\h\\l>CK^^d-i{K)\~^~-^ (26) 
for A > {2/rKY, where tk is the inradius of K. 

Proof. Recall that B'^{y, 5) C stands for a ball with center y and radius 5 > 0. We 
consider the set 



|x e i^^ : p{x, K) < e, Vo\{{B^{x, 2e) \ B\x, e)) n iT) > ^ | 



for e < rx/2. By Lemma 4 in [13], it is known that 



W^-^ {{x e dK ■ r(x) >e})> ( 1 - 



d-i 



KiVd-i{K), 



where r(x) is the radius of the largest ball that is contained in K and contains x. 
It is easy to see that x G with p{x,K) < e is in if r(proj^(x)) > e. As a 
consequence, we have 



Vol(M,) > n"-' {{x e dK : r{x) > e}) e 

> 11-^' 



(27) 



For X G Mp it holds 



|/i(a:)| > ^exp(- {a" - 2-") Xk.s") (l - exp(-2-^ Xk.s")) 



(28) 



To see (28), the underlying idea is that for every x G there is, by definition of M^, 

a set f/ C {B'^ix, 2e) \ B'^{x, e)) n K with 4(f/) = 2-'^Kde'^. Then 

¥{r]{B\x, Ae) \ f/) = 0, ri(U) > 1) = exp(-(4^ - 2-'^)\Kds'^) (l - exp(-2-'^ Xt^ds'^)) 

and for this event the effect of adding x to the point process is larger than ^^||-. 
Combining (27) and (28), we obtain 

ll/ill? = A /" f\{xfdx>\ [ f\{xfdx 



> \-,K^^d-i{K)e ^exp(-2(4'^ - 2-'^)A«:,£'^) (l - exp(2-'^ Xk^s'')) 



2^ 

and the choice e = leads to 

^2 



HAH? > ^exp(-2(4'^-2-'^)/s:,) (l - exp(-2-'^ Ka))' mVd-i{K)X-^--^ . 



□ 
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Remark 3 An inequality as (26) cannot hold for all A > as the following consider- 
ation shows: 

We fix a compact convex set Kq with & K and a compact window^W^ D K and set 
Kr = rK = {rx : x G Kq} for r > 0. We define the random variable PV(iy) as 

PV(W^) =Vo\{{y eR'^ : p{y,W) < \\y - x\\ \/x e t] nW^'}) . 

A short computation proves E PV{Wy < oo. Then, it holds 

VarPV(ir^) = KFV{Kr)^ - Vol(ir^)2 < EPV(K^)2 
< (1 -exp(-AVol(/s:,)))E PV(iy)^ 



For r — 7- 0, the right hand side has order r'^, whereas Vd-i{Kr) is only of order r 



d-l 



5 Proof of Theorem 1.1 



In this section, we use our abstract central limit theorem Theorem 3.1 to prove Theorem 



1.1 



Since it follows by (25) that the condition (11) is satisfied, it remains only to check 



(12), which requires 



VarPV(fs:)' VarPV(fs:) 
We show that for every tt G Hij 



— 7- as A — )■ oo. 



Ma = A'^I / . . . / \{fi*fi*fj* fj)Ayu • • • , y\n\)\dyi . . . dy\^\ 



converges to zero as A — ?■ oo at a sufficiently high rate such that the inequalities (16) 



and (17) in Proposition 3.4 imply that condition (12) is satisfied. 
We define functions (?„ : (R'^) " — )• M, n G N, as 

Qnixi-i . . . , Xn) = max < diam(xi, . . . , a;„), max p{xi, dK) 

I j=l,...,n 

where diam(a;i, . . . , Xn) stands for the diameter of xi, . . . , Xn- Using this notation, we 
can state the following upper bound for /„: 

Lemma 5.1 Lei Xi , . . . , x„ G M'^ and 6 = gn{xi, . . . ,Xn)- Then 

|/„(xi, . . . , xOI < 4v exp(-A/€,574'^) =: 
nlA 

Proof. Without loss of generality we can assume p{xi,dK) = 6 or p{xi,X2) = S. For 
the first case, let x = ^{xi + projg^(a;i)), where projg^(a;i) is the projection of Xi on 
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the boundary of K. If the projection is not unique (this can happen for Xi G K), it 
does not matter which of the points is taken. Then, it holds 

^ < p(z/,x) <\5< p{y,dK) for all y G B''{x^,S/4). 



Hence, it follows from Lemma [4. 1| and a straight forward computation as in the proof 
of Lemma 14.21 that 



+ - [ Fiz{y,v)^B\y,p{y,dK)\\))dy 

< I nz{y.v)iB\y,\\y-Xi\\))dy 

+ - I ¥{z{y,ri)iB\y,\\y-i\\))dy 

= 4vexp(-A/c,574'^). 
n\X 

In the case p(a;i, X2) = 5, we replace x by X2 and obtain the same bound. □ 

We prepare the application of the coarea formula by showing the following properties 
of Qn ■ 

Lemma 5.2 a) Qn is a Lip schitz function with ||V(?„|| > 1 almost everywhere, 
b ) There is a constant Q > only depending on the dimension d such that 

d-i 



VJ'''-\9-\5)) < Y,>,d-MK)5^~'dnd5''-\>,d5''r-^ (29) 

i=0 

d-l 
d-l 



1=0 



Proof. gn{xi, . . . , Xn) is always given by the distance of two points or by the distance 
of a point to the boundary of K. If we move one of these points exactly in the opposite 
direction v of the second point or the boundary of K, the directional derivative is 
VvQuixi, . . . ,Xn) = \ \v\ \ and 

|Vi,5'„(xi, ...,Xn)\ < ||V5'„(xi, . . . ||-i;|| 

implies ||V(y'n(2;i, . . . ,Xn)\ \ > 1 and thus a). 
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For the proof of b) we consider the same situations as in the proof of a). If there 
are two points Xi,Xj G M'' such that p{xi,Xj) = S, Xi must be in {dK)s, xj in a sphere 
around Xt with radius 6 and the remaining n — 2 points must be in a ball with radius 
S and center Xi. If p{xi,dK) = S, Xi must be in the set {y E : p{y,dK) = 6} 
and the remaining points are in a ball with radius 6 and center Xj. Combining these 
considerations with the Steiner formula (see Theorem 2.2.4 in |ljy) yields (29). □ 



For / = 1,2,3,4 let nii^yi, . . . ,y\Tr\) C {yi, . . . ,y\TT\} be the new variables that occur in 
the l-th function of (/j * fi* fj * fj)n and let \ni \ stand for the number of these variables. 
We set r = g\^\{yi, y\jr\) and 

^1 = 5'kii(7ri(yi, . . .,y\n\)), . . . ,^4 = g\7T4\Myi, ■ ■ ■,y\n\))- 

Since 7ii{yi, . . . ,y\^\) C {yi, . . . ,y\^\}, it is easy to see that 6i = g\^^\{ni{yi, . . . ,y\^\)) < 
g\n\{yi, • • • , y\TT\) = r for / = 1, 2, 3, 4. If there is a yj with p{yj, dK) = r, we have at least 
two li,l2 e {1,2,3,4} such that yj E tt,, (?/i, . . . , ?/|^|) and yj E ni^{yi, . . . ,yi^i), which 
implies = = r. The other case is that there are yj-^ and yj^ such that p{yji,yj2) = 
r. If there is a / G {1,2,3,4} with 2/1,1/2 G ni{yi, ■ ■ ■ ,y\Tr\), it follows directly 61 = r. 
Otherwise, vr G Hij implies that we have a yj^ and /i,/2 G {1,2,3,4} with yj^^yj^ E 
7r«i(yi, . . . ,2/|,r|) and yj.„yj.^ E -Ki.ijji, . . . ,y\^\). By the inequality r = piyj^^yj^) < 
piVji^yja) + P{yj3^yj2), it follows max{6i^,6i,} > max{p{yj^,yj^), p{yj^,yj^)} > r/2. 
Hence, it holds r/2 < max;=i ... 4 5; < r. 

Together with the coarea formula Theorem 2.2, Lemma 5.1 and Lemma 5.2, we 
obtain 



< aI"I-^ 



fiiSi)fi{52)fji5s)fj{Si)dyi ...dy\ 



^1 



(Rd)|7r| 



fi{5i)fi{52)fj{h)fj{5A)\\^9U\\\dyi ■ ■ ■ dy\^\ 



Jgu^ir) 
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i^.(5i)/^(52)/,(53)/,(54)Hl'^l^-^(%i, . . . , y\^\))dr 



poo 

/ exp(-A/s:,r78^)Q|7r| V K,_MK)iKyp-'^ dr. 
-^0 i=0 



By substitution and the definition of the Gamma function, we have 

d-l 

I Irtt^ 1 1 I 

My < 



(z!)2(j 
16crf|7r 



I " ^ roo 

^5^«:,_,V,(ir)8W'^— ^A^-^ / exp{-XKy/8''){\Ky/8''p-'^dr 



1=0 
d-l 



i^^«:,_.r.(ir)A^-^8W'^-V/ / exp( 



rflp^ E .._.F.(ir)A^-8W-^., ^r(|.| - Vrf). 



i=0 
d-l 



j=0 



Thus, each Mx has the order A d and, by Proposition 
most the same order. As a consequence of Theorem 1.2 



3.4 



all Rij and Ri have at 



VarPV(-ft') has the order 
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A ^ d , which means that 
and 




and 



Var PV(X) 
' Ri 



have a order less or equal than A 



— )• as A — )■ oo. 



VarPV(/s:)' VarPV(fs:) 
Now all assumptions of Theorem |3.1| are satisfied for the volume of the Poisson-Voronoi 



approximation and Theorem 1.1 is a direct consequence. 



Remark 4 In Theorem 1.1 and Theorem 1.2, we assume that the approximated set 
K is convex. But the convexity is not necessary for the construction of the Poisson- 
Voronoi approximation such that it is a natural question, if one can extend our results 
to more general set classes. The convexity assumption is only needed to bound the 
volume and the surface area of the parallel sets {dK)r by the Steiner formula in the 
proofs of Theorem 1.2 and Theorem 1 1 . 1 1 and to apply Lemma 4 from [TS] in the proof 
of Lemma 4.6, Hence, one can extend the results to compact sets M C M"^ that satisfy 
the following additional assumptions: 



(SI) There are constants c 



, d, depending on M such that 



Vol((9M),) < ^c^r' and 'H'^-\d{{dM)r)) <Y, 



i=l 



i=l 



with {dM)r = {x eR'^ : p{x, DM) < r} for r > 0. 



(S2) It holds 



liminf Vol(M,.)/r > 



r-s>0 



with Mr = {xeM^ : p{x,M) < r,Yo\{{B'^{x,2r) \ B'^{x,r)) f] M) > 
Assumption (SI) allows us to bound the volume and the surface area of the parallel 



sets {dM)r by a kind of Steiner formula. In the upper bound in Theorem 1.2, the 
intrinsic volumes must be replaced by the constants cj^], z = 1, . . . , rf. Our proof of the 
lower bound in Theorem 1.2 requires assumption (S2), which replaces a rolling ball 

Then the constant C in ([T| and the lower bound for A 



result for convex sets from 
depend on the limit inferior in (S2). 

Since (SI) and (S2) are obviously true for convex sets, they still hold for polyconvex 
sets and, of course, for all polytopes. 
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